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Abstract. We present accurate mid-infrared observations of the disk around the young, bona-fide brown dwarf CFHT-BD- 
Tau 4. We report GEMINI/T-ReCS measurements in the 7.9, 10.4 and 12.3 fim filters, from which we infer the presence of a 
prominent, broad silicate emission feature. The shape of the silicate feature is dominated by emission from 2 //m amorphous 
olivine grains. Such grains, being an order of magnitude larger than those in the interstellar medium, are a first proof of dust 
processing and grain growth in disks around brown dwarfs. The object's spectral energy distribution is below the prediction 
of the classical flared disk model but higher than that of the two-layer flat disk. A good match can be achieved by using an 
intermediate disk model with strongly reduced but non-zero flaring. Grain growth and dust settling processes provide a natural 
explanation for this disk geometry and we argue that such intermediate flaring might explain the observations of several other 
brown dwarf disks as well. 
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1. Introduction 

Circumstellar disks play an importa nt role in the f ormation 
and early evolution of stars (see, e.g. lShu et alJI 19871) . as well 
as in the subsequent formation of planetary systems (see, e.g. 
lLissaueifl993.) . Similarly to young stars, infrar ed excess emis- 
sion was recently identified from brown dwarfs dComeron et al 



that of the T Tauri disks? Can brown dwarf disks potentially 
form planetary systems? 

Recent studies investigated the mid-infrared photome- 
try from about 15 brown dwarf disks aiming to d etermine 



the disk struct ure (see, e . g. iNatta ^Testi 200 11: lApai etal 
et al.ll2003l: iMohantv et al...2004^ 'Walker etal 



' I2OOO'; M uench et a l. 2001': ' Liu et alJ2 003": Jav awardhana et al 
[2003a). The interpretation of this excess emission as a proof of 
substantial amount of dust confined into a disk, was strongly 
supported by millimetre-wavelength observations estimating 
the total mass of two brown dwarf disks being 0.4-6 Mjup 
llKIein et alJ2003h . 

The detection of brown dwarf disks sets at least two very 
important questions: How does the disk structure compare to 
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2004 1. The two basic models used are the flat (Adams e t al 



1988h and flared (opening angle increasing with the radius, 
Chiang & Goldreichlll997l) disks. Due to their geometry the 
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flared disks absorb a larger fraction of the stellar light. This 
additional energy is re-emitted at mid- and far-infrared wave- 
lengths leading to a far-infrared bump often observed in the 
spectral energy distribution of T Tauri and Herbig Ae disks. 

The flat and flared disk model predictions deviate strongly 
longwards of 15 fj.m. In spite of the lack of observations in this 
regime, there is mounting evidence that both disk structures can 
be found in brown dwarf disks. 

The similarity of disk structures and stellar mass/disk mass 
ratios between young stars and brown dwarfs emphasizes the 
exciting question whether planet formation can take place in 
brown dwarf disks. The major steps leading to planet for- 
mation start by grain growth and dust settling. While the 
early phases of such grain coagulatioi i has been recently ob- 
serv ed in the d isk atmosphere (e.g., 'Przvgodd a et al.l 120031: 
ivan Boekel et al..2003: Meeus et al..2003.) and in the disk mid- 
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plane (e.g.. ICalvet et alJ l2002t pTesti et alJ l2003t iNatta et alJ 
I2OO4I) of young low- and intermediate-mass stars, no such ob- 
servation could be performed for the much fainter brown dwarf 
disks. 

In this letter we focus on the disk around the spectro- 
scopicall y confirmed, non- a ccreting brown dw arf CFHT-BD- 
Tau 4 (iMartfn et alJ l200lt Ijavawardhana et aL .2003b. R. A. 
04^ 39" 47.3'* Dec. +26^ 01' 39" J2000) being the best- 
understood such o bject with rneasure ments in the near-infrared 
I Liu et ali .2003>: iMartin et al.1 l200ll) . mid-infrared ISOCAM 
I Pascucci et alJl2003l) and sub-milli metre regimes (IKlein et al.l 
12003 1). Based on these photometrv Pascu cci et alJ 12003) car- 
ried out the most detailed and systematic analysis of disk mod- 
els. They excluded all but three models as possibilities: the 
single- or double-layered flat disk or the flared disk with in- 
ner rim. These two models differ in the far-infrared and in the 
presence of the 9.7 fim sihcate feature. As a continuation of 
this effort we present new, accurate mid-infrared photometry 
for CFHT-BD-Tau 4 obtained by the recently installed T-ReCS 
camera. The filter bandpasses probe the region of the 9.7 fim 
silicate emission feature. Modeling the shape of this feature al- 
lows us to determine the disk structure and test the dust grain 
properties. 

2. Observations, Data Reduction and Results 

The observations have been carried out on Jan 2, 2004 using the 
T-ReCS mid-infrared detector mounted on the Gemini South 
8m-telescope in service mode. The plate scale was 0.09"/pixel 
with a field of view of about 29"x22". Three filters were 
used with central wavelengths (50% transmission wavelengths) 
of 7.9 fim (7.39-8.08 //m), 10.38 fim (9.87-10.89 jum) and 
12.33 yum (1 1.74-12.92 yL(m). 

In order to eliminate the high thermal background the usual 
chopping/nodding technique was used with throws of 10". The 
total on-source integration times were 12, 7 and 10 minutes for 
the filters 7.9, 10.4 and 12.3 /im. The data reduction has been 
carried out by using self-developed IDL routines. Each chop- 
ping pair has been inspected manually to ensure the exclusion 
of the frames suffering from strong detector signatures, such as 
sinusoidally modulated background noise. During this proce- 
dure we rejected about 10% of the chopping pairs. Because the 
object is undetected in the individual chopping pairs, our man- 
ual selection does not bias towards stronger or fainter source 
flux. After the manual inspection the high-quality frames have 
been averaged. 

The CFHT-BD-Tau 4 observations suffered from imperfect 
tracking, which resulted in somewhat blurred images. The rel- 
atively bright source, however, could be identified also on sub- 
groups of the images, allowing us to track and correct its po- 
sition during the exposure. This correction proved that CFHT- 
BD-Tau 4 is neither extended nor elongated on the T-ReCS im- 
ages and allowed more accurate photometry. 

The ffux calibration is based on the calibrator HD 92305 
(7 Cha), which was observed several times during the night. In 
order to derive accurate flux de nsities for the calibr ator we in- 
tegrated its spectra as given by 'Cohen et alj ill999h within the 
50% transmission levels of the T-ReCS filter set. The flux mea- 



surements have been carried out by using aperture photo me- 
try w ith the IDL-adaptation of the DAOPHOT routine ( Stetsoiil 
Il987l) . To increase the accuracy of our photometry, we tested 
different aperture annuli. We found that the fluxes in the final 
two negative and one positive beams of the given source show 
consistent, slowly changing values for aperture radii between 
0.9"- 1.35". The fluxes obtained in this aperture range with 1- 
pixel steps have been used to determine the count rates from 
the source. These count rates have been converted to absolute 
flux densities by comparing them to the mean of three count 
rate measurements of the standard star obtained by the corre- 
sponding aperture sizes. 

The mid-infrared emission of this brown dwarf disk is un- 
resolved at a spatial resolution of 0.3". Our images prove that 
no adjacent mid-infrared object could have influenced the pre- 
vious coarse-resolution ISOCAM photometry. 

We find the following flux densities for CFHT-BD-Tau 4: 
F7.9 = 42+10 mJy, Fm.4 = 51.4+2.6 mJy andFi2.3 = 47.1+2.4 
mJy. The final fluxes are the mean values over the aperture 
range, while the photometric error bars given are the maximum 
deviations from the mean value inside the aperture range com- 
bined with the calibration uncertainties. The large error bar on 
the 9 measurement reflects the low atmospheric transparency 
and the short integration times. The excess emission at 10.4 fim 
filter is undoubtly marking a prominent 9.7 fim silicate feature. 



3. IVIodels and Interpretation 

In this section we compare the new T-ReCS observations to 
model predictions aiming to explore the dust properties and 
disk geometry. First, we will use the mid-infrared data to con- 
strain the most dominant dust species. Then, we will apply this 
information to a set of self-consistent disk models with differ- 
ent geometries. 

3.1. Dominant Dust Species 

We model the dust composition by simply decomposing the ob- 
served mid-infrared (T-ReCS and ISOCAM) flux densities into 
the sum of a black body and a modified black body, represent- 
ing the continuum emission and the optically thin component 
responsible f or the emission feature . A similar approach ha s 
been used in lBouwman et al.l ( 1200 ll) and lMeeus et alJ ( l2003l) . 
and for the details of the fitting procedure we refer to those 
articles. The weighting factors of the two components and the 
common temperature have been fitted by using the Levenberg- 
Marquardt minimization. By introducing different optical con- 
stants to the optically thin component we tested dust grains of 
diverse types and sizes as possible dominant species. We note, 
that the pure black body component can also represent dust 
grains without spectral features in the mid-infrared regimes, 
such as carbonaceous materials. The weighting factor of the 
modified black body is proportional to the mass in the optically 
thin regime. 

In our model we used silicates wit h olivine compo- 
sition (FeMgSi04, Dorschn er et alJ Il995l) . the most abun- 
dant interstellar silicate species, as well as amorphous silica 
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dSpitzer & Kleinmarll960l) . For both species we fitted the ob- 
servation using spherical grains with radii between 0.1 and 
100 fim. Based on the;^^^ comparison we found, that the 2.0 yum- 
sized olivine grains provide the best fit to the observations. The 
best-fit temperature was 490 K and the total dust mass in the 
optically thin regime equaled to 0. 1 lunar mass. 

In order to put a qualitative constraint on the mass ratio 
of the large/small grains, we included a second optically thin 
component. We find, that the lower limit for the mass ratio of 
the 2.0 /im grains to the 0.1 fim ones in the optically thin disk 
regime is m2.o/?Mo.i > 6. 

3.2. Self-consistent disk models 

The fitting procedure described above has no self-consistent 
physical basis, but it provides a good indication for the dom- 
inant dust species. We confirm this i nformation bv usi ng self- 
consistent disk models described in IPullemond et al.1 fc oOl). 
which are improved versions of models by Chiang & Goldreichi 
In these models the radiation field of the central object 
together with the dust properties determine the disk flaring an- 
gle. 

In IPascucci et al.l ll2003l) we already applied these mod- 
els for CFHT-BD-Tau 4 and compared several disk geome- 
tries. We excluded all but three models: single-layer flat disk, 
two-layer flat dis k and flared disk with inner rim (m4, m5 
and m9 in Pascucc Tet alJl20o'3l) . The main parameters of those 
and the following models are: brown dwarf luminosity Lbd - 
0.1 Lo, mass Mbd = 0.075 Mq, temperature Tef j = 2800 K; 
the dust opacities ar e those of O.lyum astronomical silicates 
jPraine & Leelll984 . The outer radius was set to 50 AU, but 
our observations do not constrain this parameter. The pres- 
ence of the 9.7 /zm silicate feature arises from an optically thin 
disk layer (e.g., disk atmosphere). As the single-layer flat disk 
model lacks any optically thin component, this model can be 
immediately excluded. When comparing the m9 flared and m5 
flat disk models, we find that none of them matches well the 
new observations as long as the 0.1 jum astronomical silicates 
are used as dust species. Thus, successfully fitting the detailed 
mid-infrared data requires the identification of the dust species 
and its typical grain size in addition to the disk structure. 

Here we use the results of Sect. 13.11 namely, our guess of 
2.0 jum olivine grains as the dominant dust species. When we 
apply 2.0 jum-sized amorphous olivine grains in the disk mod- 
els, we find that the two-layer flat disk model (m5) reproduces 
the observations better than the flared disk model. The best fit, 
however, is reached by applying a disk model intermediate be- 
tween the flat and flared disk models: as discussed in the Sect.|3 
we artificially reduced the disk height by a factor of 4 to ac- 
count for dust settling effects. 

As an illustration of the influence of the main parameters 
and the stability of our results in Fig. [l] we show the best- 
fit models with flat and flared geometries in combination with 
interstellar-like 0. 1 //m astronomical silicate grains, as well as 
the intermediately flaring disk model and 2 //m grains. 

Throughout the modeling procedure we also varied the in- 
chnation the inner disk radius i?,„, the exponent a of the sur- 



[x." 40 
20 



. 0.1 astron. sil. 


' / 


^ — 


, flared 


/ flared . 
y + rim 

C 










flat 


. intermediate flaring 




2 /<m olivine 














I'/ 










2 ixm astron. sil. 




6 


8 


10 


12 14 

m] 


16 



Fig. 1. Flat, flared and intermediate disk models compared to 
the observed spectral energy distribution of CFHT-BD-Tau 4. 
The upper panel has been calculated by using 0. 1 astronomical 
silicate grains, while the lower panel applies 2.0 /vm olivine and 
astronomical silicate grains. All model parameters are identi- 
cal to those in J'ascucci et al. (2003), except for the reduced 
flaring parameter of the intermediate model. The 850 //m and 
1.3 mm observations are not shown here but fitted well, sim- 
ilarly to Pascucci et al. (2003). The best fit is achieved by the 
intermediately flaring model in combination with 2yum amor- 
phous olivine grains. 

face density power-law r", the presence of an inner rim and 
self-shadowing effects, the latter two being self-consistently 
calculated. The best solutions are found around the values of 
/ ^ 0°, Of = -1.9, Rin = 3.5 /?BD, and no inner rim. 

Independently of the grain size estimate obtained by de- 
composing the emission feature, we also tried to fit the observa- 
tions by using different dust species and geometries. However, 
no other configuration than a disk with reduced flaring and 
moderate-sized (2 /im) olivine grains could reproduce the mea- 
surements. 

4. Discussion 

4.1. Disk structure and Dust settling 

The sparsely sampled spectral energy distributions of the 
brown dwarf disks observed until now have been generally 
modeled by either flat or flared disk geometries. The flared disk 
models are based on the assumpti on of efficient mixing of ga s 
and dust everywhere in the disk ("Chiang & Goldreich"l922l)- 
In this picture the dust particles are kept above the disk mid- 
plane by the turbulent gas, which in turn is thermally coupled 
to the dust grains absorbing the incident stellar radiati on. Grain 
growth, however, reduces the dust-ga s coupling (Dubrulle et alJ 
Il995t iDullemond & Dominikll2004 and accelerates the sink- 
ing of dust particles toward the disk midplane. Thus, the dust 
disk converges towards the flat disk geometry. Reduced disk 
scale heights have been argued for in the case of several disks 
around low- and intermediate-mass stars (see, e.g. lChiang et alJ 
.2001.) . 
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Our modeling yields two major results: 1, the optically thin 
disk regime is dominated by at least intermediate-sized (~2/im) 
dust grains; 2, the dust disk's scale height is strongly reduced. 

These results are fully consistent with the evolution of 
disks through grain growth and subsequent settling, as outlined 
above. Up to now, this process has been only rarely observed in 
circumstellar disks and never in a disk of a brown dwarf. The 
multi-wavelength observations of the relatively bright disk of 
CFHT-BD-Tau 4 provide the first evidence that the similarity 
between young stars and brown dwarfs extends to circumstel- 
lar dust processing and the subsequent evolution of the disks. 

We note here, that the picture of grain growth and dust 
settling might provide a natural explanation not only for 
CFHT-BD-Tau 4 but for several ot her brow n dwarf disks 
with flat disks, such as Cha Hal dApai et ; 



^l2Q02h. GY5 
)m dNatta & Testil 



/Mohantv et al.'2004") and further 9 objects from 
boOL-Nat ta et al. 200^). 

In Sect. 13.11 we estimated the mass in the optically thin 
regime to be ~0.1 lunar mass. The same method when applied 
to the sample of Herbig Ae stars presented in Bouwman e t al] 
results in masses of ~ 100 lunar masses, about a thou- 
sand times larger than that of CFHT-BD-Tau 4. If these disks 
would be optically thin, their mid-infrared emission ratio would 
be similar to their mass ratio, i.e. ~ 50. The larger difference, 
however, is naturally explained if the emission originates from 
the atmosphere of an optically thick disk. In this case, the total 
emission scales with the surface area contributing to the 9.7 jjm 
silicate feature: The radial dust temperature distribution in the 
disk atmosphere approximately scales as T(r) cc T, x (r)"?, 
where r is the radius, r» is the stellar temperature, assuming the 
absorption coefficient /c to be oc v. Given the typical Herbig Ae 
stellar temperatures being ~4 times larger than the correspond- 
ing brown dwarf temperatures, any given temperature around a 
brown dwarf is reached at ~ 30 times smaller radii compared to 
the Herbig Ae stars. This translates to a surface ratio of ~ 10^, 
the same order of magnitude as the fitted mass ratios. 
Our best-fit model 



similarly to tho se of Natta & Testil 
(l200lh . iMohantv et all tOOA IWalker et al. L2004,) and others 
— require an inner disk radius of Rin - 3.2 Rbd- This value is 
similar to t he typic al inner disk truncation radii deduced for T 
Tauri stars dShu et al., 1994) and might be defined through dust 
sublimation at temperatures of ~ 1500K. 

4.2. Grain growth in a substellar disk 

Our observations demonstrate that the disk atmosphere of 
CFHT-BD-Tau 4 is dominated by 2.0 fim grains, an order of 
magnitude larger than the average grains commonly found in 
the interstellar matter. Further modeling indicates that the grain 
growth is accompanied by dust settling. These two processes 
are th e initial, essential steps of planet formation (LissaueJ 
Il993l) . 

Studies of coeval T Tauri disk systems showed that the 
dust evolutionary stage is only a weak functio n of time (see, 
e.g. |Przvgoddaet'aDl2003t iMeeus et alJl2003l) . On the other 
hand, the rate of dust processing is likely to provide a use- 
ful way to parameterize the evolution of an accretion disk 



through protoplanetary disk into a debris disk. Comparing the 
derived ot o n/fftn i ratio to those e stimated for disks of Herbig 
Ae stars dBouwman et al.l I2OOII) provides some insight into 
the stage of the dust evolution in CFHT-BD-Tau 4. The val- 
ues range from m2.o/mo,i < 0.04 for the unprocessed dust 
in the Galactic Centre up to > 57 for the evolved disk of 
HD 100546. The ma.o/'wo.i > 6 ratio of CFHT-BD-Tau 4 indi- 
cates that this brown dwarf di sk is close in dust e volution phase 
to the 5 Myr-old HD 104237 fcradv et alJ2004l) . The fact, that 
comets Halley and Hale-Bopp show the presence of smaller 
grains (ma.o/'Wo.i < 3), might suggest that they formed in a dust 
evolutionary phase - but not necessarily in time - prior to that 
of CFHT-BD-Tau 4. 

We directly derive the feature stren gth from the observa- 
tions following IPrzveodda et alJ (l2003h and compare its ~ 1.6 
value with the T Tauri s tar samples of*Mee us et al ] il2003h and 
przygodda et aO (l2003l) . This model-independent comparison 
again shows CFHT-BD-Tau 4 among the disks with the most 
processed dust. 

The young age of CFHT-BD-Tau 4 (~1 Mvr. iMartm et alJ 

l200lh stands intuitively in contrast to the evolved stage of its 
dust disk. There are three factors which may resolve this ap- 
parent contradiction: First, large differences in dust processing 
timescales have been spotted among T Tauri and Herbig Ae 
stars, indicating that dust processing can proceed at very differ- 
ent rates. Second, the over-luminosity of CFHT-BD-Tau 4 may 
be explained by an unresolved binary brown dwarf. Dynamical 
perturbations could then influence the dust coagulation and 
settling timescal es. Hint s for this behavior have been found in 
T Tauri dis ks (Meeus et al. 2003) and other brown dwarf disks 
dSterzik et al.a20 04). Third, if it is a binary system, then the age 
of CFHT-BD-Tau 4 is likely to be underestimated. 

The only previous detection of silicate emission from a 
brown dwarf disk was found for the p Ophiucus object GY310 
(Mohanty et al. 2004). In that case the shape of the feature hints 
on small dust grains, but the data quality do not allow any firm 
conclusion. Spectroscopic surveys with the new Spitzer Space 
Telescope will be able to derive the accurate dust composi- 
tion for a large number of brown dwarf disks revealing how 
frequent grain growth processes are. However, to follow the 
grain growth to scales larger than a few micron, high-resolution 
(sub)millimetre observations are necessary, possibly by using 
interferometric facilities such as the Sub-millimeter Array or 
ALMA. 



5. Summary 

The main conclusions of this work are the following: 

1 . Using the T-ReCS/Gemini instrument we detected thermal 
infrared emission at 7.9 jjm, 10.4 fim and 12.3 fim from the 
disk of CFHT-BD-Tau 4. 

2. The disk displays a prominent silicate emission peak which 
proves the existence of an optically thin disk layer 

3. We compare the shape of the silicate emission feature to 
laboratory spectra emerging from different dust species and 
grain sizes. Using two independent methods we find that 
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the emission feature is dominated by 2 /urn amorphous sili- 
cate (olivine) grains. 

4. Simple semi-analytical models explain the observations by 
a two-layered flared disk with a reduced scale height. This 
transitionary disk model between flat and flared geometries 
is consistent with dust settling models. 

5. The dominance of 2 fim grains and the reduced disk flaring 
prove that grain growth and settling occurred in the disk 
of the ~1 Myr-old CFHT-BD-Tau 4. Based on the derived 
mass ratio for the 2 fim to 0.1 yum grains the stage of dust 
processing in this brown dwarf disk is similar to that of the 
well-known Herbig Ae star HD 104237. 

In view of these results we conclude, that the dust grain coag- 
ulation and dust settling, the initial steps of planet formation, 
take place in the disk of the brown dwarf CFHT-BD-Tau 4. 
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